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I 

Abstract 

_  The  parametric  instability  growth  rate  of  ion  acoustic  and  Langmuir 

waves,  driven  unstable  by  two  uniform  pumps  near  the  Bohm-Gross  frequency,  is 
calculated  as  a  function  of  pump  amplitudes  and  frequencies.  TVo  instability 
mechanisms  can  be  identified:  one  corresponds  to  the  usual,  single  pump 
parametric  instabilities  (decay  and  oscillating  two  stream)  while  the  other 
is  similar  to  that  found  in  the  Mathieu  equation.  The  interaction  between 
these  two  mechanisms  results  in  a  non-monotonic  dependence  of  the  growth  rate 
on  the  pump  amplitudes  and  frequencies:  both  cancellation  and  enhancement  are 
obtained  for  various  values  of  the  parameters.  An  analytic  study  of  the 
relevant  dispersion  relation  using  Hill's  method  is  complemented  by  numerical 
studies  in  both  the  frequency  and  time  domains. 
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U  INTRODUCTION 

The  growth  rates  of  Langmuir  and  ion  acoustic  waves  driven  unstable  by  a 
high  frequency  uniform  electric  field  E0(t)  =  Ejcosu^t  +  E2cosa)2t  in  a  uniform 
unmagnetized  plasma  are  calculated  for  "pump"  frequencies  cdj  and  w2  close  to 
the  Bohm-Gross  frequency  ook  =  (o>p2  +  3Tek2/m)1/2,  with  the  difference  fre¬ 
quency  A  =  -  o>2  being  of  order  of  the  ion  acoustic  frequency  £lk  =  kcs. 

While  earlier  work  on  this  problem1  considered  the  threshold  for  instability, 
i.e.  the  minimum  value  of  (Ej2  +  E2^)  required  to  amplify  a  noise  signal,  we 

examine  here  the  growth  rate  as  a  function  of  ui],  u>2,  Ej  and  E2.  Experimental 
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observations  on  both  experimental  and  ionospheric  plasmas  have  shown  a  com¬ 
plicated  dependence  of  the  wave  amplitudes  on  pump  frequencies.  Although  most 
steady-state  experiments  sample  the  nonlinearly  saturated  state  of  the  insta¬ 
bility,  it  is  important  to  determine  the  properties  of  the  early  linear  growth 
stage.  The  present  work  addresses  this  question  and,  in  addition  to  providing 
the  necessary  conceptual  framework  for  a  future  nonlinear  saturation  theory, 
predicts  various  nontrivial  features  which  may  be  useful  in  understanding  the 
experimental  observations.  Since  most  experiments  operate  well  above  the 
threshold  levels,  the  damping  of  both  Langmuir  and  Ion-acoustic  waves  is 
neglected  here;  it  then  suffices  to  use  the  warm  fluid  theory  rather  than  a 
kinetic  theory  model.  The  calculations  can  be  extended  to  include  damping 
terms  but  at  the  expense  of  introducing  more  parameters. 

The  characteristics  of  the  usual  single  pump  parametric  excitation  are 
well  known.  When  the  frequency  matching  condition  wi  -  <dr  “  flk,  is  satisfied 
there  is  a  "decay  instability",  consisting  of  the  decay  of  the  pump  wave  at  u>j 
into  a  Langmuir  wave  and  an  ion  acoustic  wave,  both  having  wavenumber  k.  For 
any  finite  pump  amplitude  Ej  there  is  a  range  of  u>]  around  wk  +  Jlk  for  which 
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the  growth  rate  y  of  the  daughter  waves  is  positive.  Equivalently,  for  given 
<D1  in  that  range,  there  is  a  threshold  for  Ej.  Of  course,  when  the  frequency 
matching  condition  (Oj  =  <0^  +  is  satisfied  exactly,  the  threshold  drops  to 
zero,  when  damping  is  neglected.  When  to ^  =  u^,  the  so-called  oscillating 
two-stream  instability  (OTSI)  occurs,  but  here  we  shall  discuss  primarily  the 
decay  instability.  All  of  these  properties  are  immediate  consequences  of  the 
fourth  order  differential  equation  whi  'h  governs  the  time  evolution  of  the 
spatial  Fourier  transform  of  the  ion  density  n^(k,  t).  Since  the  coefficients 
of  this  equation  are  constants,  its  Laplace  transform  yields  a  simple  disper¬ 
sion  equation,  a  quadratic  in  u)2,  whose  solution  gives  the  results  stated, 
previously . 

With  two  pumps,  the  coefficients  in  the  differential  equation  for  n^  are 

not  constant.  Instead,  they  become  periodic  functions  of  t,  with  frequency 

A  =u»]_  —  o)2,  resulting  in  a  differential  equation  which  resembles  the  well- 

known  Mathieu  equation,  albeit  of  higher  order.  Physically,  this  Mathieu-like 

character  arises  from  the  ponderomotive  force  at  the  beat  frequency  A  driving 

the  ion  acoustic  waves.  As  might  be  expected,  the  solutions  of  this  equation 
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have  properties  analogous  to  those  of  the  Mathieu  equation.  Since  the  Laplace 
transform  of  this  equation  leads  to  an  infinite  set  of  coupled  equations  for 
the  quantities  n^(k,  <d  ±  £A),  where  l  is  an  integer,  the  dispersion  relation 
takes  the  form  of  the  vanishing  of  an  infinite  determinant,  an  equation  which 
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we  solve  for  w  using  Hill's  method. 

The  behavior  of  the  resulting  solutions  can  be  described  in  terms  of  two 
separate  instability  mechanisms,  one  similar  to  the  usual  single  pump  paramet¬ 
ric  instability,  the  other  analogous  to  that  found  in  the  Mathieu  equation. 
Specifically,  if  <*>j  is  near  +  %,  if  Ej  exceeds  the  single  pump  threshold 
Es  and  if  E2  is  of  order  Ej ,  we  recover  the  usual  decay  instability,  except 
for  certain  values  of  o>2,  where  the  growth  rate  y  vanishes  or  is  slightly 


enhanced.  This  is  Illustrated  In  Fig.  4  where  the  real  and  Imaginary  parts 
of  a)  for  the  growing  wave  are  plotted  as  a  function  of  <1)2  for  fixed  u>j  and 
Ej  =  E2.  In  the  figures,  and  In  subsequent  sections,  we  use  the  dimensionless 
quantities  va  =  and  ga  =*  (eEa/Te)(kkp)~1/2(M/16m)1/lf  to  characterize 

the  pump  frequencies  and  amplitudes.  Fig.  4  illustrates  some  general  features 
of  the  equal  pump  amplitude  case  with  Ei  =  E2  =  E  >  Es.  Although  y  s  Im  u>  is 
equal  to  the  single  pump  growth  rate  for  most  values  of  u>2,  it  vanishes  when 
vi  +  v2  =  0  or  vi  +  v2  =  2  and  it  is  enhanced  slightly  when  v2  =  Vj  -  2/N, 
where  N  is  an  integer,  or  when  v2  a  0  (corresponding  to  the  OTSI  driven  by  the 
second  pump).  As  we  explain  in  more  detail  in  Sec.  IV,  these  results  are 
representative  of  cases  where  the  parametric  decay  instability  tends  to 
dominate  the  behavior  but  is  modified  by  the  Mathieu-like  effects.  (Further 
details  concerning  the  results  shown  in  Fig.  4  and  in  the  other  figures 
mentioned  in  this  introductory  section  are  given  in  Sec.  IV.). 

If  oil  is  near,  but  not  exactly  equal  to  +  ill  (vj  =  1)  and  Ej  is  below 
the  single  pump  threshold  Es  corresponding  to  this  value  of  uj,  then  only  the 
Mathteu-ltke  effects  can  produce  instability,  as  illustrated  in  Fig.  5a,  5b, 
and  5c  where  E  *  Ej  =  E2  is  successively  increased,  but  remains  below  Es.  In 
Fig.  5d,  where  E  exceeds  Es,  one  sees  a  combination  of  the  two  effects:  the 
single  pump  decay  mechanism  dominates  for  v2  far  from  -1,  while  the 
Mathieu-like  effects  give  enhanced  growth  near  v2  *  -1  and  zero  growth  for 
finite  intervals  of  v2  above  and  below  -1. 

For  unequal  amplitudes,  there  is  a  complicated  interplay  between  the 
Mathieu  and  decay  mechanisms,  the  resulting  behavior  depending  on  the  ratio 
E2/Ej  and  also  the  ratio  Ej/Eg.  If  Ej  <  Eg,  the  decay  instability  does  not 
occur  (unless  v2  =  1)  but  the  Mathieu  instability  appears,  as  illustrated  in 
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Fig.  6a  ,  the  growth  rate  at  first  increasing  as  E2  increases  (Fig.  6b),  and 
eventually  decreasing  (Fig.  6c).  If  Ej  is  well  above  Es>  the  second  pump  may 
simply  modify  the  decay  instability  growth  rate,  as  shown  in  Fig.  7a  (details 
in  Sec.  IV)  for  E2  =  E^/2,  or  it  may,  for  larger  E2,  actually  suppress  the 
growth  rate  entirely  over  a  finite  interval  of  V2  between  the  Mathieu  and 
decay  instability  regions  as  in  Fig.  7b. 

If,  instead  of  fixing  uij  and  varying  a>2  we  keep  u>2  constant  and  vary  u>j, 
we  observe  a  similar  interaction  of  the  two  instability  mechanisms  as  shown  in 
Fig.  8.  In  general  a  mixture  of  the  two  mechanisms  is  most  likely  to  occur 
when  vj  =  1,  maximizing  the  growth  rate  of  the  decay  instability  and, 
simultaneously  V2  -  -1  so  that  Vj  -  Vj  =  v  =  2  which  corresponds  to  the 
strongest  Mathieu-like  instability.  Most  of  our  attention  has  been  focused  on 
this  "mixed  regime". 

Since  the  parameter  space  (Ej,  E2,  u>2)  is  four  dimensional,  surveying 

it  is  greatly  facilitated  by  having  an  approximate  solution  of  the  dispersion 
equation.  Judicious  truncation  of  the  infinite  determinant  yields  a  simple 
approximate  dispersion  equation  (a  biquadratic  in  id)  which  gives  close 
agreement  with  the  exact  results  and  also  provides  a  simple  means  of 
understanding  the  propeties  of  the  solution  of  the  exact  dispersion  relation 
displayed  in  Figs.  5  through  8.  In  addition,  this  approximate  dispersion 
relation  can  be  used  to  determine  the  boundaries  of  the  stable  and  unstable 
regions  in  the  (Ej,  E2,  <V(,  0*2 )  space  as  illustrated  in  Fig.  3  for 
two-dimensional  cross  sections  (E2  vs.  V2  for  fixed  Ej,  Vj;  Ej  vs.  Vj  for 
fixed  E2,  ^2*  ®*2  vs*  V1  f°r  and  v2)*  Although  this  approximation  can 

be  strictly  justified  only  when  V£  *  -V2  =  1,  It  proves  to  be  valid,  in  fact, 
over  a  fairly  broad  range  of  parameter  space,  as  illustrated  in  Fig.  9 
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As  an  alternative  to  the  Laplace  transform  approach,  we  have  also  solved 
the  differential  equations  in  the  time  domain,  using  numerical  integration,  a 
procedure  which  avoids  the  various  approximations  used  in  solving  the 
dispersion  relation.  Taking  the  Fourier  transform  of  the  solutions  for  nj(t) 
and  ne(t)  then  gives  directly  the  "line  shapes"  n^co)  and  ne(u>)  which  would  be 
observed  experimentally  if  the  nonlinear  saturation  mechanism  were  independent 
of  frequency.  Although  these  quantities  can  also  be  calculated  from  the 
solutions  of  the  Laplace  transformed  problem,  that  approach  would  give 
spectral  peaks  corresponding  to  all  roots  of  the  dispersion  equation,  growing, 
decaying,  or  stable,  whereas  in  the  time  domain  calculation  (and  in  the 
experimental  situation)  the  growing  waves  dominate.  The  time  domain  solutions 
also  show  clearly  the  modulational  effects  which  result  from  the  occurrence  of 
two  unstable  roots  of  the  dispersion  equation. 

The  theoretical  model  used  is  presented  in  Sec.  II,  together  with  a 
derivation  of  the  differential  equation  for  nj(k,  t)  and  a  discussion  of  its 
similarity  to  the  Mathieu  equation.  An  exact  solution  of  the  frequency  domain 
equations  is  given  in  Sec. Ill,  where  we  also  show  that  a  judicious  truncation 
of  the  infinite  determinant  leads  to  simple  expressions  for  the  frequencies 
and  growth  rates  and  for  the  boundaries  of  the  stable  and  unstable  regions  in 
parameter  space.  The  results  of  these  numerical  calculations  and  a 
discussion  of  the  various  features  are  presented  in  Sec.  IV.  Section  V  gives 
the  solutions  in  the  time  domain  for  both  growing  and  stable  waves. 

Conclusions  are  presented  in  Sec.  VI. 
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II.  MODEL  EQUATIONS  AND  DISPERSION  RELATION 

We  consider  a  uniform,  unmagnetized  ion-electron  plasma  with  a  uniform 
''pump"  electric  field: 

E0(t)  =  costojt  +  E2  cosu^t  (1) 

Pump  depletion  is  neglected,  so  _Ej  and  Ej?  are  constant;  and  are  assumed 
parallel;  and  the  pump  frequencies  <Uj(j  =  1,  2)  are  near  the  electron  plasma 
frequency  u>p.  Since  we  are  interested  in  growth  rates  well  above  threshold, 
wave  damping  is  neglected. 

The  fluid  equations,  for  each  species,  01  =  e,i, 

+  7  *  (n  v  )  =  0  <2> 

3t  a-a 

3v 

~a  +  v  •  7va  +  7pa/nama  =  qaE/nio,  (3) 

3t  0 

are  linearized  about  the  oscillating  motion  due  to  Id 

na  =  no  +  nla»  Pa  =  poa  +  Pla*  Xxx  =  2ioa  +  Zla»  Zoa  =  claSo/,ma 


This  gives 

nla  *  ^^}Loa* ^ la  ^ioa*^^nla  =  ^  Pia^ma  —  QL  ~ 

if  terms  of  order  v£a  are  neglected. 

Using  Poisson's  equation;  neglecting  the  zeroth  order  ion  velocity  (vcj) 
since  the  ions  do  not  respond  to  the  high  frequency  field;  setting  T^  =  0 
(i.e.,  taking  the  limit  of  large  Te/T t);  Fourier  analyzing  in  space;  and  sep- 
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arating  the  electron  density  Into  high  frequency  (to  of  order  <op)  and  low  fre¬ 


quency  parts,  ne  =  neh  +  ne£,  we  obtain 

••  2  • 
neh  +  “k  neh  =  "ik  <vo  ne*> 

2 


Ope  (nejj,  -  nj)  +  k2Te  nei/m  =  ik  (v0  n^)* 


(■>) 

(6) 


**  2  2 
ni  +  “pi  ni  =  “pi  ne£ 


(7) 


Here  we  have  used  pe  =  aneT  with  a  =  3  for  the  high  frequency  and  a  =  1  for 
the  low  frequency  equation.  Also,  we  set  me  =  m,  it^  =  M,  u>£  =  top2  +  3k2Te/m 
and  neglect  terms  of  order  m/M.  Inserting  (6)  into  (7)  and  approximating 

“pi  11  “  “pe/(“pe  +  k2Te/M^I  =  k2Te/M  =  k?  cs  =  we  8et 

"i  +  nk  ni  =  -(ike/M)(E0neh)£  (8) 

4h  +  “k  neh  =  (ike/m) (Eonj)  (9) 


Note  that  the  right  hand  side  of  (8)  corresponds  to  the  usual  ponderomotive 
force.  Finally,  a  modulational  representation 


neh  =  f+(0  exp(-itokt)  +  f_(t)  exp(itokt)  (10) 

where  f+  are  slowly  varying  functions,  lf+/f+l  <<  tok,  gives  the 
equations: 

32n 

_ L  +  nA  =  -  (ix/2) ( f _( t )A( t )  +  f+(T)A*(T>]  (11) 

3x2 

8f+  ,  -  ( xM/ 1 6m ) 1 / 2  niA  (12) 

3t 

3f-  -  (xM/16m) niA* 

3  x 


(13) 
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where 

2  2 
A  =  l  Aj  exp  (-iVjT),  Xj  =  (ekEj/mu>p  ), 

vj  =  (u)j-a»k)/flk,  X  =  (kD/k)2,  and  t  =  Jlkt. 

The  remainder  of  the  paper  is  concerned  with  the  properties  of  the  solutions 
of  (11)  through  (13).  Of  course,  these  equations  are  only  valid  for  small  Xj 
since  we  have  dropped  terms  of  order  v02. 

Before  discussing  the  solution  of  (11)  through  (13),  we  note  that  solving 
(12)  and  (13)  for  f+  and  substituting  the  result  into  th  o.  deromotive  force 
expression  on  the  right  side  of  (10)  gives  terms  of  the  ..oral  form  A  A*  n^. 

(The  actual  analysis,  given  later  in  this  section,  acto  ,  *nvolves  differen- 

j!  tiating  (10)  to  obtain  a  sixth  order  equation  for  n^;  the  approximate  discus- 

# 

sion  in  this  paragraph  is  only  meant  to  illuminate  the  physics  involved.)  If 
j  there  is  only  a  single  pump,  A  A*  is  constant  so  the  equation  for  nj  has  con¬ 

stant  coefficients;  the  only  effect  is  to  change  the  eigenfrequency  from 
w  =  1  to  a  new  value  which,  for  pump  amplitudes  above  the  decay  or  OTSI 
thresholds,  becomes  complex.  However,  with  two  pumps,  A  A*  contains  also 
oscillating  terms  of  frequency  v  =  vj  -  V2  and  it  is  this  oscillating  pondero- 
motlve  force  term  in  the  ion  density  equation  of  motion  which  is  responsible 
for  the  new  effects  arising  with  two  pumps. 

A  direct  method  of  solving  these  equations  is  to  use  the  Laplace  trans¬ 
form 

00 

n^(w)  =  /  dr  nj(T)exp(lu)r) 

o 

(where  w  is  measured  in  units  of  Uk)  which  leads  to  the  set  of  coupled  equa- 


t  ions 
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with 


Y_(w)ni(<»)-v )  +  X(w)ni((u)  +  Y+fuOn^iD+v)  =  I(to) 
2 

X(o>)  =  0)2-1  +  |  g^2v.  /( 0)2-v^2) 

j  =  l 


Y  (u>)  =  (vl+v2)gjg2 
±  2(oi±Vj)(o>+V2) 


(14) 


(15) 


2  1/2  2 
gj  =  X  (Mx/16m)  Aj 


Eq.  (14)  also  follows  directly  from  Eq.  (25)  of  Arnush  et  al.  1 ,  derived  from 
the  Vlasov  equation,  provided  we  take  the  fluid  limit  for  the  e(k_,  o>)  in  their 
equation.  The  right  hand  side  of  (14)  involves  the  initial  conditions  on 
n i ( t )  and  its  derivatives,  but  we  can  simply  set  it  equal  to  zero  in  finding 
the  dispersion  relation.  Thus,  vanishing  of  the  determinant  D  of  the  coupled 
equat Ions  ( 1 4) 


D  =  det 


(-1) 

(-1) 

,(°) 

(0) 

(0) 

- 

X 

(1) 

(1) 

(1) 

Y+ 


=  0  (16) 


yields  the  dispersion  relation,  where  x(°)  =  X(u>+nu>)  and  similarly  for  Y+(n). 
(Note  that  the  matrix  of  eq.  (16)  is  tridiagonal.) 

Before  discussing  the  solution  of  (16)  we  examine  the  Mathieu-like  equa¬ 
tion  derived  from  the  system  (11)  through  (13).  Differentiating  (11)  twice 

and  (12)  and  (13)  once  allows  us  to  eliminate  f+  and  obtain  a  fourth  order 
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equation  for  n{ 

(4)  (1)  (3)  (2)  (1)  (1) 

S  n^  -  S  n^  +  ( S+H )  n£  —  S  +  (CS2+H)  nj  —  0  (17) 

1/2 

where  C  =  (Mx/m)  (x/8),  the  upper  script  (j)  indicates  the  jth  derivative 
with  respect  to  t  and  S  and  H  are  periodic  functions  of  vt  =  (v>i-V2)t: 

S  =  2[  [  Vj  A ^ 2  +  a1A2(vi-V2)cosvt) 

H  =  Vj3  Aj2  +  A1A2V|V2(Vi+V2/cosVTJ 

Because  of  the  S2  term  in  the  coefficient  of  n^,  the  Laplace  transform  of 
this  equation  couples  n^(u>)  not  only  to  n^(u>±v)  but  also  to  nj(a>±2v),  i.e. 
we  get  a  5  term  recursion  relation  rather  than  (14).  For  this  reason,  it  is 
more  convenient  to  work  instead  with  a  sixth  order  equation  for  nj  whose 
coefficients  involve  only  cosvx  and  sinvr  but  not  the  harmonics  of  v.  This 
equation  can  be  obtained  by  operating  on  (11)  with  the  two  operators 


Lj  =  (d/di)2  +  Vj2  j  =  1,  2  (18) 

and  using  Eqs.  (12)  and  (13)  to  simplify  the  terms  on  the  right  side.  The 
result  is 


[LiL2L  +  (g12vlL2  +  g22v2L|)j  nt  = 

-(Vi+w2)  (L2-  Lj+{exp(ivi)n^}  +  L2+  L1“{exp(-ivx)n1} 


(19) 


where 

L  =  (d/dt)2  +  1  Lj  *  =  (d/dr  ±  iVj) 


(20) 


Of  course,  the  Laplace  transform  of  Eq .  (19)  gives  just  Eq .  (14). 

Equation  (19)  is  an  interesting  generalizaton  of  the  Mathieu  equation, 
which  we  can  write  as 


1 1 


d-^n  +  +  g2n  cos vt  =0  (21) 

3t2 


Indeed,  the  Laplace  transform  of 
(14)  but  with  X,  Y+  replaced  by 
X(a»)  =  u)2  -  l 

Y  (w)  =  -  g2/ 2 


(21)  gives  a  system  of  equations  identical  to 


(22) 


and  an  analysis  of  that  system  leads  to  the  usual  Mathieu  stability  diagram. 

A  comparison  of  (19)  and  (21)  is  instructive.  If  the  right  side  of  (21),  i.e. 
the  part  with  periodic  coefficients,  vanishes,  the  dispersion  equation  reduces 
to  X(u>)  =  0,  which  has  only  the  stable  solutions  w  =  ±1.  On  the  other  hand, 
even  if  the  right  hand  side  of  Kq.  (19),  i.e.,  the  part  with  periodic  coeffi¬ 
cients,  vanishes,  the  resulting  dispersion  equation  (14),  X(w)  =  0,  has  both 
stable  and  unstable  roots.  Indeed,  with  g2  =  0  (which  makes  the  right  side  of 
(19)  vanish)  X(o>)  =  0  is  just  the  usual  dispersion  for  single  pump  parametric 
instabilities 


(m2-l)  (o)2-v}2)  +  g  j  2v  j  =  0  (23) 

This  has  unstable  roots  for 

812  >  ( v i2“l )2/4v i  or  0  <  -vt  <  gl2  (24) 

corresponding  to  the  parametric  decay  and  the  OTSI,  respectively.  More  gener¬ 
ally,  X(w)  =  0  has  unstable  solutions  for  given  vj,  V2  if  exceeds  certain 
threshold  levels.  This  corresponds  to  one  mechanism  for  instability,  which 
may  be  considered  as  a  straightforward  extension  of  the  usual  parametric  in- 
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stabilities  (decay  and  OTSI),  and  is  quite  different  from  the  situation  in  the 
ordinary  Mathieu  equation,  where  there  is  no  instability  in  absence  of  the 
term  with  a  periodic  coefficient. 

The  second  (Mathieu-lik.e)  mechanism  is  evident  when  gj  is  below  the 
threshold  for  instabilities  arising  from  X(u)  =  0,  since  Eq .  (19)  for  g2  *  0 
can  still  have  instabilities  due  to  the  terms  on  the  right  hand  side  with 
periodic  coefficients.  This  instability  mechanism  is  clearly  analogous  to 
that  associated  with  the  Mathieu  equation,  where  instabilities  arise  only  from 
the  g^  term  in  (21).  In  general,  for  arbitrary  gj  and  g2  we  have  the  presence 
of  both  instability  mechanisms. 
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I  It.  ANALYSIS  OF  THE  DISPERSION  EQUATION 

In  this  section  we  solve  the  dispersion  equation  (15)  using  Hill's  meth¬ 
od.  We  also  discuss  an  approximation  which  provides  some  physical  insight  and 
which  proves  to  be  quite  accurate  in  the  parameter  regime  of  greatest  inter¬ 
est.  The  approximation  is  based  on  an  expansion  in  the  quantity  g  =  (gig2)*^ 
which  we  treat  as  a  small  parameter. 

For  arbitrary  g  we  define  a  new  determinant  D  obtained  by  dividing  D  by 
its  diagonal  elements 

00  (n) 

D(u>)  =  D(oi)  n  X  (w) 

n=-» 


D(u>)  = 


where 

(n)  (n)  (n) 

W  =  Y  (u>)/X  (u>) 


is  of  order  g2.  if  g  is  small,  we  can  expand  D  in  powers  of  g: 


00  (n)  (n+1) 

D<«)  =1+1  W+  (m)  W_  (w)  +  0(g)  (26) 

In  the  limit  g  *>  0,  D  =  1  but  D  can  still  vanish,  namely  if  w  is  a  root  of 

(n) 

X  (w).  For  later  use  we  define  the  quantities  rf,  i  =  1  to  6,  as  the  roots 


of  X(w),  for  arbitrary  gf  and  g2 
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X(rt)  =  r i2  -  1  +  1  gj2Vj(r12-Vj2)-1  =  0  (27) 

j 

In  general,  the  roots  r^  have  no  simple  physical  significance,  but  we  note  the 
following  properties.  If  gj  *  0  and  g2  =  g  =  0,  then  Eq.  (27)  Is  a  biquadrat¬ 
ic  In  oj  whose  roots  give  the  usual  single  pump  parametric  instabilities  as 
illustrated  in  Fig.  1  for  gj  =  .346.  (The  points  A,  P,  M  shown  in  Fig.  1 
are  used  later  in  the  discussion  of  Fig.  5.).  In  the  stable  region 
| v  j  —  1 1  >  g^  there  are  two  real  roots,  the  ion  acoustic  mode  with  a>  =  1, 
and  what  we  may  call  the  pump  idler  mode,  with  to  =  Vj.  If  gj  and  #2  are  both 
non-zero,  Eq.(27)  is  a  cubic  in  ui2  whose  roots  are,  for  small  gj,  close  to 
u>2  =  1  (normal  ion  acoustic  modes,  present  even  if  gj  0)  and  to  u>2  =  ooj2 
and  u)2  =  u)22  (pump  modes  or  idlers,  which  have  no  physical  significance  in  the 
limit  gj  *  0). 

We  now  consider  the  roots  of  D(u>)  for  small  g.  Since  the  second  term  In 
6  is  of  order  g^,  6  can  vanish  only  if  one  or  more  of  the  X^n^(to)  is  small, 
i.e.  if  (a)  +  nv)  is  near  l,  corresponding  to  an  ion  acoustic  resonance.  Of 
particular  interest  is  the  case  when  two  of  the  vanish  simultaneously, 

which  can  happen,  with  v  *  0,  if,  for  two  integers  nj,  and  n2» 

ui  f  njv  =  1  a)  +  02^  =  -1  (28) 

i.e.  if 

v  *  2/N  (29) 

where  N  is  an  integer.  Of  special  interest  is  the  case  of  "double  resonance" 
where  one  of  the  terms  in  D  involves  the  product  of  two  large  W  factors,  e.g. 
when 

a)  =  -1  <i)  +  v  -  1  (30) 

which  requires  v  =  2,  N  =  1.  The  most  interesting  results  are  obtained  in 

this  case,  which  can  be  understood  on  physical  grounds  as  follows.  From  the 
differential  equation  (19)  we  see  that  the  ions  are  driven  by  a  ponderomotive 
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force  of  frequency  v,  so  a  low  frequency  wave  at  to  gives  rise  to  another,  at  fre¬ 
quency  u-v.  Both  waves  can  be  normal  modes  if  to  =  1,  v  =  2  so  that  to  -v  =  -1, 
giving  the  condition  of  "double  resonance".  For  N  >  1  this  double  reso¬ 
nance  cannot  occur  so  the  effects  of  the  second  pump  are  less  pronounced. 

However,  the  case  N  =  2,  v  =  1  is  of  some  interest.  A  wave  at  to  =  1  then 
gives  rise,  through  the  ponderomotive  force,  to  a  wave  at  to  =  1-v  =  0,  and 
while  this  is  not  a  normal  mode,  it  does  correspond  to  an  OTSI  mode.  The 
thresholds  for  this  case  are  discussed  in  Arnush  et  al.1 

Having  seen  where  the  most  interesting  effects  are  likely  to  occur,  we 
consider  the  exact  solution  of  the  dispersion  equation  (15),  i.e.  without 
assuming  g  small.  The  determinant  5(to)  defined  by  Eq.  (25)  has  the  following 
properties: 

a)  D(to  +  v)  *  D(u>) 

(n) 

b)  Lim  D(oj)  =  1  since  lim  W  (to)  =  0 

(i)-*oo  U»00  — 

c)  D(to)  =  D(-to)  =  6*  ( U>* ) 


(n)  (n)_ 

Since  X  (to)  vanishes  at  to  =  r^  =  r^  +  nv  where  rj  are  the  6  roots  of  X(u>), 
(n) 

D(u>)  has  poles  at  r^  .  These  will  be  simple  poles  provided  we  avoid  the 
special  v  values  where  rj  -  r^  =  pv  for  some  integer  p.  Then  the  function 


K(to)  i  D(to)  -  l 
i  =  l 


“  (n)  -1 

2  b1>n(u)-ri  ) 

n--co 


(31) 


(n) 

where  b^>n  is  the  residue  of  I)  at  r^  ,  will  be  analytic  in  the  whole  to  plane 

and  since  K  +  1  as  u  *  «,  we  have  K  =*  1.  The  periodicity  of  B  implies  that 
(n)  (m) 

the  residues  at  r^  and  r^  must  be  the  same,  i.e.  that  =  bj  indepen¬ 


dent  of  n.  Finally, 
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“  (n)  "I  -  -1 

l  )  =  l  v-1  [n-(o)-ri)/v]  =  (Tt/v)cot[it(u-ri)/v]  (32) 

n=-°°  n=-°» 

SO 

6 

D(u>)  =  1  +  1  (bjn/v)  cot[TT(oi-rj[)/v]  (33) 

i=l 

The  symmetry  of  both  5  and  X  under  ui  -*■  -w  means  that  If  we  arrange  the  roots 
r^  so  that  r^+3  =  -r j ,  i  =  1,2,3,  then  bj+3  =  -bj  and  the  dispersion  relation 
takes  the  form 

3 

f)(u>)  =  1  +  1  (b^  ir/v){cot(ir(tu~ri)/v]  -  cot [ w(a)+r^ ) /v ] }  =  0  (34) 

i  =  l 


Since  (34)  gives  the  dependence  of  0  on  w  in  explicit  form,  it  is  easy  to  de¬ 
termine  the  roots  of  D  once  the  bj  are  known  as  functions  of  v.  Note  that  the 
method  fails  when  v  =  0  since  the  single  poles  r^11)  converge  into  a  single 
point  giving  rise  to  an  essential  singularity.  Therefore  the  neighborhood  of 
v  =  0  is  excluded  in  our  numerical  calculations.  However,  v  =  0  implies  that 
both  pumps  have  the  same  frequency  and  the  corresponding  growth  rate  is  expec¬ 
ted  to  be  that  of  a  single  coherent  pump  whose  amplitude  is  the  sum  of  the  two 
pump  amplitudes. 

The  r^'s  are  poles  of  D,  so  the  residues,  b^,  are  given  by 

bj  =  lim  D(u))(a)-ri)  (33) 

w+ri 

These  infinite  determinants  can  conveniently  be  evaluated  using  iteration: 

®n+l  _  a  _  an,n+lan+l ,n 

Dn  n+1  ,n+l  (Dn/D^j)  (36) 

where  Dn  is  the  approximate  value  of  6  obtained  from  an  n  x  n  truncation  and 
an  +  i,n  +  l»  an,n  +  l»  an+l,n  are  t*ie  elements  of  the  (n+l)th  row  and  column.  Al- 
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though  some  of  the  an  m  are  singular  at  u>  =  rj  due  to  the  vanishing  of  X(u>) 
factors  which  appear  in  the  denominator,  (u)-r^)all  m(io)  is  always  Iloilo. 

The  iteration  scheme  converges  fairly  well;  typically,  for  values  of  v  bigger 
than  0.5,  10  iterations  suffice  to  give  an  accuracy  of  0.1%. 

The  roots  of  D  obtained  by  this  procedure  are  plotted  as  functions  of  the 
parameters  vj,  V2,  etc.  as  discussed,  with  illustrative  examples,  in  Sec. 

IV.  To  understand  the  considerable  structure  which  results,  the  following 
approximate  treatment  of  the  double  resonance  case  proves  helpful.  For 
v  =  2,  u)  -  1  we  have,  to  order  g\  keeping  only  the  largest  resonant  terms, 

D  =  1  +  W+(u>-v)  W_(ui)  =  1  +  Y  ( u>-v)  Y_(u))/X(o>-v)  X(w)  =  0  (37) 

This  equation,  which  is  equivalent  to  approximating  D  by  a  2  x  2  determinant, 
gives  an  eighth  degree  polynomial  in  a),  when  rat ionalized,  but  it  can  be  re¬ 
duced  to  a  biquadratic  in  the  limit  vj  =  -V2  a  1.  Since 

X(o>)  =  (m2-i)  +  g12/2(o)-v1)  -  g22/2((D+v2) 

X(w-v)  =  {(id-v)2  -  1}  -  gj2/2(u+v2)  +  g22/ 2(iu-v j ) 

Y+(io-v)  =  Y_(to)  =  (v  j+v2)(g  jg2)/2( u)-Vj  )  (u»+v2)  , 
the  change  of  variable 

u>  =*  v/2  +  y  (38) 

gives  a  cubic  in  y2, 

(y2-d2)(y2-a2)(y2-b2)+(g 12+g22)b(y2+ad)-(g12-g22)vy2-( l/4)(g12— g22)2  =  0 


(39) 
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where  a  =  (<5i  -  62)/2,  b  =  (6j  +  <S2)/2,  6|  =  V]  -  1  and  62  =  v2  +  *  are  small 
quantities  in  the  double  resonance  regime  and  d  =  v/2  +  I,  with  d  =  2  in  that 
regime.  If  the  growth  rate  y  and  the  shift  in  the  normal  mode  frequency, 
(Reu>-1),  are  much  smaller  than  1,  as  is  the  case  for  small  amplitude  pumps, 
then  |  y  1 2  <<  ^2  for  v  ~  2  and  equation  (39)  becomes 


y4-2Ay2  +  B  =  0 

where 

A  =  -i  (a2+b2+(g 12+g22)b/d2-(g 12~R22) v/d2) 
B  =  a2b2  +  (gi^-g2^)2/4d2  -  (g i2+g2^)ab/d 


(40) 


(41) 


From  (38)  and  (40)  we  have 

oi  =  v/2  +  [A  +  (A2-B)1/2]1/2  (42) 

as  a  convenient  closed  form  approximation  for  the  roots  of  D(oj)  which  is  valid 
for  Vj  =  -V2  -  1  and  small  values  of  gj.  (Note  that  (42)  gives  four  roots 
since  each  square  root  can  have  either  a  positive  or  negative  sign.) 

The  rather  complicated  dependence  of  the  growth  rate  on  the  pump  fre¬ 
quencies  and  amplitudes  as  determined  from  the  numerical  solution  of  (34)  can 
be  understood  in  a  fairly  simple  way  from  an  examination  of  (41)  and  (42), 
which  actually  provide  a  good  approximation  even  when  6j  and  62  are  not  small. 
From  (42)  we  can  see  that  there  are  two  disjoint  conditions  for  instability 


a)  B  <  0 

b)  B  >  A2 
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(A3) 
(44) 

and  the  boundaries  between  the  stable  and  unstable  regions  of  parameter  space 
are  then  determined  by  the  locii  B  =  0  and  B  =  A2,  as  shown  in  Fig.  2.  We 
first  note  that  the  number  of  unstable  roots  of  (42)  varies  from  0  to  2 
according  to  the  signs  of  A,  B  and  A2  -  B.  The  various  possibilities  are  sum¬ 
marized  in  Fig.  2. 

We  may  say  that  region  a)  represents  the  generalized  Mathieu  instability 
since  Re  u>  is  locked  to  one  half  of  the  (ponderomotive  force)  driving  fre¬ 
quency,  v/2,  just  as  in  the  lowest  unstable  mode  of  the  usual  Mathieu  equa¬ 
tion.  Similarly,  we  may  consider  region  b)  as  the  generalization  of  the 
single  pump  decay  instability  since  in  the  limit  g2  ♦  0,  V2  -1,  we  have 
B  »  A2  and  (42)  reduces  to  the  usual  single  pump  expression 

m  =  v/2  ±  A1/2  ♦  {1  +  vj  ±  [ ( 1 —v ! ) 2  -  gi2v1]1/2}/2,  (45) 

1/2 

unstable  if  g  >  gs  =  (l-Vj)Vj- 

The  mapping  of  the  stability  boundaries  on  the  physical  space  of  the 
parameters  can  be  obtained  by  examining  the  surfaces  B  =  0  and  B  =  A2  which 
characterize  the  different  regions  of  Fig.  2.  Since  there  are  four  indepen¬ 
dent  parameters  (vj,  \>2,  g\  and  g2)  it  is  convenient  to  fix  two  of  these  and 
plot  the  curves  B  =  0  and  B  =  A2  in  the  plane  of  the  remaining  two  parameters. 
Examples  of  such  plots  are  shown  in  Fig.  3  with  the  respective  regions  of 


Fig.  2  identified. 


In  this  section  we  analyze  the  roots  of  (34)  as  functions  of  the  parame¬ 
ters  Vj,  \>2»  and  ^2*  Numerical  calculations  of  both  growth  rates  and  fre¬ 


quencies,  are  presented  and  the  approximate  form  of  the  dispersion  equation 
(42)  is  used  to  analyze  the  resulting  structure. 

The  new  effects  resulting  form  the  presence  of  the  second  pump  are  most 
evident  at  resonance,  i.e.  ,  when  one  of  the  pump  frequencies  differs  from  the 
Bohm-Gross  frequency  by  approximately  the  ion  acoustic  frequency.  We 
shall  therefore  consider  two  cases: 

i)  =  1.  v2  arbitrary 

ii)  \>2  -  “1*  arbitrary 


i)  vi  =  1 

We  begin  with  the  case  of  equal  amplitude  pumps  gj  =  g2  =  g»  with  g  well 
above  the  singLe  pump  decay  instability  threshold,  gs  =  (vj-1 )vj~l/2  =  6j.  Of 
the  various  roots  of  (34)  we  select  the  ion  acoustic  wave,  i.e.,  the  one  with 
Y  =  Im  oj  >  0  which  has  Re  <v  =  1  (in  units  of  fi^).  In  Fig.  4  we  show  y  and 
Re  w  as  a  function  of  V2  for  vj  =  1 . 1  and  gj  =  g?  =  R  =  0.346.  (The  region 

\>2  ~  vi  is  excluded  from  the  plot  for  the  reason  stated  in  Sec.  III.  Note 
that  the  growth  rate  Increases  as  we  approach  V2=Vj,  consistent  with 
our  expectation  that  y*2ys  at  that  point). 

1/2 

For  most  values  of  V2,  we  see  that  y  -  ys  -  [vi(g2-g|)]  '  ,  the  single 
pump  decay  Instability  growth  rate,  i.e.,  the  second  pump  has  little  effect. 
The  enhancement  of  y  when  V2  =  0  is  not  unexpected,  since  the  second  pump 
could  then  produce  the  oscillating  two  stream  instabiltiy  (OTSI)  even  in  the 
absence  of  the  first  pump.  Actually,  close  inspection  of  Fig.  4  shows  that  y 
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is  enhanced  for  both  positive  and  negative  values  of  V2»  whereas  the  usual 
OTSI  instability  arises  only  for  a  pump  frequency  below  the  Bohm-Gross  fre¬ 
quency.  However,  the  most  striking  feature  of  Fig.  4  is  the  total  supression 
of  the  instability  for  \>2  =  -Vj  and  V2  =  -(2-Vj),  arising  from  the  interaction 
between  the  two  mechanisms,  parametric  and  Mathieu-like,  discussed  in  Sec.  I. 

An  understanding  of  the  structure  of  Fig.  4  can  be  obtained  from  the  ap¬ 
proximate  solution  of  (34)  given  by  (42).  For  equal  pump  powers  not  too  far 

above  threshold,  A  is  positive,  since  6j,  62  and  g[  =  g2  =  g  are  all  small 

2  2 

quantities  of  the  same  order  and  hence  the  (a  +  b  )  term  in  A  dominates. 
Therefore,  condition  (43)  for  the  generalized  Mathieu  instability  becomes, 
with  d  -  2, 

(a b)2  -  g2ab  >  0  (46) 


which  is  equivalent  to 

2  2  2 
0  <  «!  -  62  <  4g 


(47) 


In  Fig.  4,  g  >  61  so  the  right  half  of 
condition  | 6 2 1  <  | 6 j |  just  corresponds 
1  and  2  in  Fig.  4,  In  this  region,  we 
is  locked  to  v/2.  At  the  ends  of  that 

v2  “  -  V1  or  b  =  0 

vi  -  v2  =  2  or  a  =  0 


(47)  is  automatically  satisfied.  The 
to  the  region  between  the  nulls  marked 
find  that,  as  expected  from  (42),  Re  <0 
interval  we  have  =  ±62,  l*e* » 

(48) 


Since  (48)  is  independent  of  the  pump  amplitude,  these  are  stable 
points  of  the  system  and  can  not  be  excited  even  if  the  pump  amplitude  is  in¬ 
creased  (subject,  of  course,  to  the  small  pump  amplitude  assumption  which  un- 
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derlies  our  whole  analysis). 

We  note  that  the  nuL Ls  at  1  and  2  in  Fig.  4  correspond  to  a  certain 
symmetry  in  the  frequency  spectrum.  For  V2  =  ~vl«  the  pumps  are  symmetrical¬ 
ly  placed  above  and  below  the  Bohm-Gross  frequency,  resulting  in  a  cancella¬ 
tion  similar  to  that  which  occurs  if  a  harmonic  oscillator  is  driven  by  equal 
amplitude  pumps  symmetrically  located  above  and  below  its  resonant  frequency. 
The  existence  of  this  null  has  also  been  noted  by  Fejer  et  al.^  For 
V1  -  v2  =  2,  corresponding  to  the  null  at  2,  the  symmetry  manifests  itself 
at  low  frequencies,  as  follows.  A  low  frequency  fluctuation,  at  frequency  u), 
beating  with  the  pump  at  w j  gives  rise  to  a  sideband  at  -  u),  the  interac¬ 
tion  being  strongest  when  the  sideband  is  resonant,  i.e.  ,  when  u)j  -  oi  =  or 
to  =  v j .  The  sideband  at  toj  -  <o,  beating  with  the  second  pump,  produces  a  low 
frequency  oscillation  at  (oij  -<o)-<02  =  v-  a>=v-Vj  =  -v2*  Rouble  reso¬ 
nance  occurs  when  both  of  the  low  frequency  signals  are  near  the  ion  acoustic 
frequency,  Sl|<,  i.e.  =  -V2  =  1.  When  -  V2  =  2,  we  have  vj  -1  =  1  -  (-V2), 
i.e.,  the  two  low  frequency  signals  are  located  symmetrically  above  and  below 
ilk  ai'd  there  is  a  cancellation.  Ibis  same  situation  arises  when  both  pumps 
are  above  the  Bohm-Gross  frequency,  vi  -  v2  ~  1  ar|d  vl  +  v2  =  2  or 
vj  -  1  =  1  -  V2.  The  resulting  cancellation  accounts  for  the  null  denoted  as 
3  in  Fig.  4. 

While  the  nulls  in  Fig.  4  arise  from  these  symmetries,  the  first  two  be¬ 
ing  associated  with  the  double  resonance  condition  N  =  1,  v  =  2,  the  slight 
enhancements  in  y  correspond  to  other  values  of  N.  Those  for  N  =  2  and  N  =  -1 
are  clearly  visible  (at  v2  =  V1  "  1  and  v2  ~  vl  +  2»  respectively)  and  the 
N  =  3  peak  (at  V2  -  vj-2/3)  is  barely  visible  on  the  scale  used  for  display. 
Higher  order  interactions  (N  >  3  and  N  <  -1)  would  appear  for  larger  values  of 
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g.  As  a  final  comment  on  this  equal  amplitude  case,  we  note  that  if  the  first 
pump  is  exactly  on  resonance,  vj  *  1,  then  a  =  -b;  the  two  conditions  (48)  are 
the  same;  the  two  nulls  at  1  and  2  coalesce;  and  there  is  no  Mathieu-like  be¬ 
havior. 

In  discussing  other  values  of  the  parameters,  we  shall  concentrate  on  t':e 
N  =  1  (double  resonance)  case,  where  the  most  striking  effects  occur.  Also, 
we  will  plot  Re  id  and  y  for  all  of  the  four  modes  with  Re  w  of  order  1.  In 
general,  we  plot  only  the  positive  y  values  but,  of  course,  for  each  root  of 
the  dispersion  equation  with  y  >  0  there  is  another  with  imaginary  part  equal 
to  ~y.  If  we  keep  the  pump  amplitudes  equal,  gi  =  gv  =  8  but  put  g  below  the 
parametric  instability  threshold  gs  for  the  first  pump  acting  alone,  we  obtain 
the  results  shown  in  Figs.  5a,  5b  and  5c.  In  Fig.  5a,  g  <  6j/2  and  it  follows 
from  (47)  that  the  range  of  instability  62  is  not  the  whole  interval 
~6j  <  <52  <  <5j  as  in  Fig. 4  but  instead  only  the  portion 

(<5l2  -  4gV/2  <  |62|  <  (49) 

This  corresponds  to  the  two  growth  regions  in  Fig.  5a.  In  Figs.  5b  and  5c, 

6 | / 2  <  g  <  6j.  In  this  case,  the  instability  region  for  62  expands  to  the 
whole  interval  (-6},  6j).  Finally,  in  Fig.  5d,  g  >  6j,  we  obtain  the  Mathieu 
instability  on  the  interval  (-6j,  5j)  and  the  decay  instability  when  v2  is  far 
from  -1.  The  behavior  of  both  y  and  Re  u  in  these  figures  is  consistent  with 
the  characterization  of  the  roots  of  (40)  given  in  Fig.  2.  In  Fig.  5a  as  we 
move  from  left  to  right  we  are  first  in  region  IV  of  Fig.  2,  then  in  regions 
l,  IV,  l  and  IV.  In  Figs.  5b  and  5c  we  have  that  the  central  feature  corres¬ 
ponds  to  region  I  and  the  sides  to  region  IV.  in  Fig.  5d,  again  moving  from 
left  to  right  we  go  through  regions  III,  IV,  I,  IV  and  III. 


► 
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An  alternative  view  of  the  structure  of  Fig.  I  is  as  follows.  For  Figs. 

5a,  5h  and  5c  where  g[  <  gs  (for  a  given  value  of  V|)  there  would  he,  in 
absence  of  the  second  pump,  two  stable  modes,  the  ion  acoustic,  with  tu  =  1 
(point  A  in  Fig.  1)  and  the  "pump"  mode,  with  w  =  Vj  (point  P  in  Fig.  1).  The 
ponderomot ive  force  with  frequency  v  may  interact  with  either  of  these,  insta¬ 
bility  occuring  if  the  beat  wave  resulting  from  this  interaction  is  also 
resonant.  This  occurs,  for  example,  when  v  -  Vj  =  vj,  i.e.  for  V2  =  -Vi,  cor¬ 
responding  to  tfie  left  hand  bump  in  Fig.  5a.  Similarly,  v  -  1  =  1  or 
v?  ~  V1  -  2  corresponds  to  the  right  hand  bump.  As  g  increases,  point  M  in 
Fig.  1  moves  closer  to  Vj  (since  the  width  of  the  instability  region  is  pro¬ 
portional  to  g)  and  the  ion  acoustic  and  pump  mode  can  be  coupled  through  the 
action  of  the  ponderomot ive  force,  resulting  in  the  y  variation  shown  in 
Figs.  5b  and  5c.  Finally,  for  g  >  gs  point  M  has  moved  to  the  right  of  Vj, 
giving  growth  for  most  values  of  V2  save  for  the  interval  where  the  Mathieu 
mechanism  dominates  and  is  stabilizing  (Fig.  5d). 

For  unequal  pump  amplitudes,  we  find  similar  dependences  of  the  roots  of 
(34).  If  the  amplitude  of  the  first  pianp  is  below  the  single  pump  instability 
threshold,  gj  <  6j,  we  obtain  the  results  shown  in  Fig.  6  for 
gy  =  .075.  For  g2  =  gj/3  (Fig.  6a),  there  are  two  disjoint  regions  of  growth 
within  the  interval  1 6 2 1  <  |  <$  [  |  whereas  for  larger  g2  these  merge  into  a 
single  region  (Figs.  6b  and  6c).  To  understand  this  structure,  we  note  that 
since  gj  <  6j,  we  have  A  >  0  and  the  condition  for  the  Mathieu  instability  is 
just  (43),  which  can  be  written  in  the  form 

6 \  “(gi  +  %2^2  <  ^2  <  512  ~  (R 1  “  R2>2  (50) 

If  g2  <  61  -  gj,  (50)  predicts  two  disjoint  Instability  regions,  as  in 
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Fig.  6a.  If  <$1  -  gj  <  g2  <  <5 1  +  gj,  then  the  left  half  of  (50)  is  automati¬ 
cally  satisfied  and  there  is  only  one  region  of  instability,  whose  boundaries 
are  within  the  interval  1 6 2 1  <  |$l|,  as  Figs*  6b  and  6c.  If  g2  >  61  +  gj, 
then  (50)  cannot  be  satisfied  and  there  is  no  Mathieu  instability,  i.e.,  we 
have  crossed  from  region  I  in  Fig.  2  to  region  IV.  The  qualitative  view  given 
previously  for  Fig.  5  also  applies  to  the  conditions  of  Fig.  6. 

If,  instead,  the  first  pump  amplitude  is  above  the  single  pump  decay 
threshold,  gj  >  6  j ,  then  A  may  be  either  positive  or  negative.  In  Fig.  7a, 
where  g(  =  1.7  6^  and  g2  =  J5 1  / 2  we  pass  from  region  I  of  Fig.  2,  when  V2  is 
near  -1,  through  region  II  (A  <  0)  and  eventually  to  region  III,  at  either 
side  of  \»2*“1.  Although  the  decay  instability  growth  rate  is  somewhat  modified 
by  the  presence  of  the  second  pump,  g2/g[  is  so  small  that  there  is  no  region 
of  V2  “here  the  instability  is  completely  suppressed.  If  g2  >  gj,  then  it 
follows  directly  from  (41)  that  A  >  0  and  instability  can  occur  only  in  re¬ 
gions  I  and  III  of  Fig.  2.  This  situation  is  illustrated  in  Fig.  7b,  where 
gj  =  1.7  6j  as  in  Fig.  7a  but  g2  =  1*1  gj.  For  these  values,  (50)  predicts  a 
single  instability  region  around  V2  =  -1,  as  seen  in  Fig.  7b,  since 
g9  <  <5i  +  gj.  As  V2  decreases,  we  pass  from  region  I  through  the  stable  re¬ 
gion  IV  and  eventually  come  to  the  decay  instability  (region  III)  near  V2  =  -1.5 
(We  pass  through  these  same  regions  as  V2  increases  from  -1  to  -0.5). 

We  see  from  (50)  that  as  the  second  pump  amplitude  increases,  we  eventu¬ 
ally  supress  the  Mathieu  instability  since  for  g?  >  gi  +  <5 1  the  condition  B  <  0 
cannot  be  satisfied.  Thus,  as  g2  increases  from  0  to  gj  +  6j,  the  maximum 
growth  rate  for  the  Mathieu  instability  (which  occurs  at  62  3  0)  for  given  Vj 
and  g 1 ,  increases,  reaches  an  optimum,  and  then  decreases,  as  shown  in  Fig.  6. 
ii)  v?  =  -1 

We  consider  here  only  the  regime  V]  >  0,  i.e.,  we  study  the  modification 
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of  t he  decay  instability  due  to  the  first  pump.  Some  typical  results  are 
plotted  in  Kig.  8.  For  comparison  the  single  pump  case  (g£  =  0)  is  shown  in 
Fig.  1.  In  Figs.  8a  and  8b  the  pump  amplitudes  are  equal,  gj  =  g2  =  g  =  .346 
so  A  >  0  and  instability  can  arise  only  from  regions  I  and  III  of  Fig.  2.  For 
V2  =  -1,  we  have 

B  =  <S  i 2  ( 6  j2  -  4g2)/16  (51) 

hence  the  Mathieu  instability  occurs  over  the  interval  1 6 j |  <  2g,  as  shown  in 
Fig.  8a.  The  vanishing  of  the  growth  rate  at  Vj  =  1  is  a  consequence  of  the 
symmetrical  location  of  the  pumps  above  and  below  the  Bohm-Gross  frequency,  as 
discussed  earlier.  When  V2  is  displaced  from  -1,  e.g.  ,  V2  =  -1.2  as  in 
Fig.  8b,  we  have 

B  =  ( 6 j  -  62  )(fil2  ~  “  4g2)/16  (52) 

Thus  the  Mathieu  instability  (associated  with  B  <  0)  occurs  for 

6!2  -  2g2  <  622  <  6t2  (53) 

This  corresponds  to  the  right-hand  and  left-hand  bumps  in  y  in  Fig.  8b.  Out¬ 
side  of  the  interval  (53)  the  Mathieu  instability  does  not  occur,  but  between 

2 

the  two  bumps  B  becomes  larger  than  A  and  we  encounter  the  decay  instability 
(region  III  of  Fig.  2)  as  evidenced  by  the  central  hump  in  Fig.  8b. 

For  different  pump  amplitudes,  we  observe  various  combinations  of  the 
Mathieu  and  decay  instabilities.  Fig.  8c  shows  a  case  where  A  >  0  (since 
g2  >  gj)  and  we  move  from  region  IV,  for  Vj  =  0,  to  regions  I,  IV,  III,  IV,  I 
and  IV  as  V[  increases  up  to  2.  For  g2  <  8|>  as  in  Fig.  8d,  we  cannot  predict 
the  sign  of  A  from  simple  arguments.  In  general  it  will  depend  on  the  values 


of  gj,  gT,  and  V[,  but  for  the  parameters  of  Fig.  Hd,  it  Is  rlt'ar,  from  the 
behavior  of  both  y  and  Re  u>,  that  as  increases  from  0  to  2  we  pass,  succes¬ 
sively,  through  regions  IV,  I,  II,  III,  II,  I  and  IV. 

From  Fig.  3  it  is  easy  to  follow  the  path  through  the  stability  plane 
corresponding  to  the  curves  of  Figs.  5  through  8.  For  example,  fig.  8d  cor¬ 
responds  to  the  dotted  horizontal  line  at  g2  =  gj/2  shown  in  Fig.  3d. 

In  all  of  the  discussions  of  this  section,  we  have  used  the  approximate 
solution  (42)  of  the  dispersion  equation  to  explain  the  results  obtained  nu¬ 
merically  from  the  exact  equation  (34).  This  is  justified  by  the  close  agree¬ 
ment  of  the  exact  and  approximate  solutions  when  and  62  are  not  to°  large. 
This  agreement  is  illustrated  for  typical  values  of  the  parameters  in  Fig.  9, 
where  the  solid  line  corresponds  to  the  the  solutions  of  the  exact  dispersion 
(34)  and  the  dotted  line  to  the  approximation  (42).  Note  that  in  this  figure 
we  have  plotted  the  growth  rates  for  all  roots,  i.e.  ,  those  with  y  <  0  are  also 
included. 
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V.  TIMK  DOMAIN  BEHAVIOR 

Although  the  stability  properties  of  our  basic  equations  (10)  through 
(12)  are  fully  described  by  the  frequency  domain  analysis  presented  in  the 
previous  two  sections,  additional  physical  insight  can  be  obtained  by  examin¬ 
ing  the  time  domain  solutions  of  these  equations.  We  use  a  fourth  order 
Runge-Kutta  technique6  to  solve  the  fourth  order  system  (10)  through  (12), 
taking  as  initial  condition  a  standing  ion  acoustic  wave  [f+(0)  =  f_(0)  =  0; 
n^(0)  =  1;  n^(0)  =0].  The  integration  is  done  using  a  step  size  At  =  0.4 
(in  units  of  *)  for  256  steps.  The  result  is  unchanged  when  At  is  taken  to 
be  0.1.  We  also  calculated  frequency  spectra  from  these  solutions  using  a  fast 
Fourier  transform  with  256  sample  points.  Insofar  as  the  nonlinear  mechan¬ 
isms  responsible  for  saturation  of  the  instability  are  weakly  dependent  on 
frequency,  these  spectra  are  representative  of  what  might  be  seen  in  an  actual 
experimental  measurement  of  the  frequency  spectrum  of  the  ion  or  electron  den¬ 
sity  fluctuations.  Of  course,  the  fine  details  of  the  spectrum,  e.g.,  the  ra¬ 
tios  of  the  various  spectral  peaks  depend  somewhat  upon  the  precise  initial 
conditions,  such  as  the  ratio  of  left-going  and  right-going  ion  acoustic  waves 
and/or  Langmuir  waves. 

Fig.  10  shows  the  results  for  a  choice  of  parameters  (vj  =  l.l,  \>j  ~  “0.8 

81  =  82  "  *346)  corresponding  to  a  point  just  to  the  right  of  null  2  in 

Fig.  4.  Since  this  is  in  region  III  of  Fig.  2,  the  dispersion  relation  pre¬ 
dicts  two  modes  with  equal  growth  rate  and  different  Re  u>.  TTie  beating  of 

these  two  modes  causes  the  modulation  in  n^(t)  shown  in  Fig.  10a.  (The  exist¬ 

ence  of  the  two  modes  is  also  apparent  from  the  plot  of  |nj(aj)|2  in  Fig,  10c.) 


The  exponential  growth  of  the  two  modes  (following  an  initial  transient  per¬ 
iod)  is  reflected  in  the  plot  of  log|ni|  in  Fig.  10b. 
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In  Fig.  II  we  show  nj(t)  for  parameters  chosen  to  correspond  to  the  nulls 
1,  2  and  3  in  Fig.  4.  In  each  case  there  is  a  modulation  but,  as  expected,  no 
net  grcwth.  When  vj  =  -V2  (as  in  Fig.  11a)  we  note  that  (11)  through  (13) 
reduce  to 


3(f  +  +  f _)/3t  =  0 


(54) 


32nj73x2  +  n^  =  -i(x/2) [f+(0)  +  f_(0)]A  cos  Vjt 

from  which  it  is  obvious  that  n^  is  a  superposition  of  oscillations  at  w  =  1 

and  id  =  Vj,  resulting  in  the  modulation  seen  in  Fig.  11a. 

In  experimental  observations,  especially  in  the  case  of  the  ionosphere, 

the  most  accessible  quantity  is  the  electron  density  ne(w)  at  high  frequencies 

(of  order  of  the  Bohm-Gross  frequency),  measured,  for  example,  by  Thomson 

2 

scattering.  In  Fig.  12  we  show  the  spectral  distributions  |ne(oj)|  obtained 
by  fast  Fourier  transform  of  the  direct  solutions  of  (10)  through  (12)  for 
equal  amplitude  pumps  gl  =  g2  =  *346  with  Vj  =  1.1  and  various  values  of  V2. 
The  locations  of  the  two  pumps  are  indicated  with  dotted  lines  and  the  ampli- 
tudes  |ne(w)|  are  arbitrarily  normalized  to  the  largest  value  found  in  this 
set,  which  occurs  in  Fig.  12g.  Of  course,  the  various  peaks  of  the  high  fre¬ 
quency  electron  density  spectrum  just  correspond  to  peaks  in  the  low  frequency 
ion  density  spectrum,  which  is  shown,  for  the  same  parameters,  in  Fig.  13; 
Indeed  it  follows  from  (12)  that 


-iu)f+(u))  =  f+(t=0)  - 


1/2 

(Mx/m)  [ A j  n1(w-v1) 


+  A2  ni(co-V2)  ]/4 


(55) 


For  vj  »  1  and  V2  far  from  -1,  as  in  Figs.  12a,  13a,  the  spectrum  is  dominated 
by  the  single  pump  decay  associated  with  the  first  pump.  As  v 2  approaches  -1 
the  spectrum  becomes  modified  (Figs.  I2b,c  and  13b, c).  For  V2  =  -1,  we  are 
in  region  I  of  Fig.  2,  so  there  is  one  unstable  root  for  Re  w  >  0  (and  another 


for  Re  u)  <  0)  in  n|(ai)  and  hence  four  (two  coinciding  at  w  =  oj^)  for  ne(a>). 
For  \>i  =  -0.88  (Figs.  I2d,  13d)  we  are  in  region  III  of  Fig.  2  and  the  two  un 
stable  modes  (for  Re  id  >  0)  in  n£(w)  result  in  a  splitting  of  the  ne(<u)  peaks 
As  \>2  continues  to  increase  (Figs.  12e  through  1 2h  and  13e  through  13h)  the 
spectrum  of  ne(oj)  again  resembles  the  single  pump  decay,  with  various  modifi¬ 
cations  in  the  line  shape. 
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VI.  CONCLUSIONS  AND  DISCUSSION  OF  RESULTS 

Our  study  of  the  effect  of  a  second  pump  on  the  parametric  instability 
growth  rate  of  ion  acoustic  and  Langmuir  waves  is  valid  for  long  wavelength 
pumps  in  a  uniform  medium,  i.e.  for  Ln  >  Xp  >  Aw  where  Ln  is  the  density  gra¬ 
dient  scale  length,  Ap  is  the  pump  wavelength  and  Aw  is  the  wavelength  of  the 
waves  excited  by  the  parametric  process. 

New  effects  appear  in  the  case  of  two  pumps  because,  in  addition  to  the 
constant  term  produced  by  a  single  pump,  which  results  in  the  usual  parametric 
instabilities,  the  ponderomotive  force  contains  also  an  oscillating  term  of 
frequency  A  =  u>j  -  012.  This  gives  the  ion  density  equation  a  character  simi¬ 
lar  to  the  Mathieu  equation,  although  the  differential  equation  arising  here 
is  of  higher  order.  We  find  that  the  interaction  between  the  Mathieu  type  of 
Instability  and  the  usual  parametric  decay  instability  (we  have  concentrated 
here  on  the  decay  instability  but  similar  results  hold  also  for  the  OTSl)  is 
strongest  when  the  condition  for  the  decay  instability  (u>j  =  1^  +  Sl^  or 
v i  *  1)  and  the  condition  for  the  Mathieu  instability  (A  =  20^  or  v  =  2),  are 
simultaneously  satisfied,  i.e.,  when  =  -V2  -  1.  The  interaction  may  be 
either  constructive  or  destructive.  Constructive  interference  is  exemplified 
by  the  fact  that  even  if  both  pump  amplitudes  Ef  and  E2  are  below  the  thresh¬ 
old  Es  for  the  single  pump  decay  instability  (in  fact,  even  if 
2  2  2 

E 1  +  E2  <  Es  )  the  Mathieu  mechanism  can  still  lead  to  instability,  as  il¬ 

lustrated  in  Fig.  5a,  5b  and  5c.  The  destructive  aspect  is  illustrated  by  the 
occurrence  of  nulls  in  y  as  a  function  of  vj  or  V2,  nulls  which  may  occur  even 
when  E |  >  E8,  as  Illustrated  in  Fig.  4  and  Fig.  7b,  and  which  may  extend  over 
a  finite  Interval  of  Vj  or  V2,  as  in  Fig.  7b.  As  discussed  in  Sec.  IV,  the  sta¬ 
bilizing  or  destabilizing  effect  of  the  second  pump  does  not  vary  monotonically 
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with  E2,  as  illustrated  in  Fig.  6  for  Kj  below  the  single  pump  threshold:  the 
growth  rate  initially  increases  with  E2,  reaches  a  maximum  and  then  decreases, 
eventually  vanishing  for  sufficiently  large  E2.  If,  as  is  sometimes  the  case, 
one  wishes  to  eliminate  parametric  instabilities,  it  appears  that  multiple 
pumps  could  be  used  if  their  parameters  are  appropriately  chosen.  Note  that 
the  suppression  of  the  instability  observed  here  occurs  with  coherent  (fixed 

phase)  pumps  and  hence  differs  from  the  use  of  broadband,  randomly  phased 

8 

pumps  . 

It  is  also  important  to  note  the  strong  dependence  of  the  growth  rate  on 
the  parameters.  For  example  Figs.  7a  and  9c  correspond  to  the  same  value  of 
the  parameters  except  for  a  slight  difference  in  the  value  of  Vj.  It  is  ob¬ 
served  that  in  Fig.  7a  the  immediate  neighborhood  of  V2  =  -1  corresponds  to 
region  I  of  Fig.  2  (i.e.,  only  one  complex  root  with  y  >  0) ,  followed  by  re¬ 
gion  III  after  a  brief  transition  to  region  2  as  [62!  increases,  whereas  the 
loop  in  y  around  \>2  =  -1  seen  in  Fig.  9c  corresponds  to  region  II  followed  by 
region  III  as  | 6 2 1  increases. 

The  solutions  of  the  dispersion  equation  discussed  in  section  IV  were  ob¬ 
tained  by  solving  the  infinite  determinant  using  Hill's  method.  However,  in 
the  parameter  regime  where  interesting  effects  appear,  an  excellent  approxima¬ 
tion  is  obtained  by  using  a  2  x  2  truncation  of  that  determinant. 

Direct  solution  of  the  problem  in  the  time  domain  is  used  to 
corroborate  the  frequency  domain  results  and  shows  the  modulational  effects 
associated  with  multiple  roots  of  the  dispersion  equation.  Moreover,  the 
Fourier  transform  of  the  time  domain  solutions  indicates  the  spectral  line 
shapes  to  be  expected  for  the  ion  and  electron  densities,  at  low  and  high  fre¬ 
quencies,  respectively.  So  far  as  the  spectral  lines  are  concerned,  the  lar- 
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gest  effect  of  the  second  pump  is  to  diminish  their  maximum  amplitudes  as  we 
approach  the  nulLs  of  Y*  but  it  also  results  in  some  fine  structure,  e.g. 
splitting  of  the  lines,  appearance  of  satellites,  etc. 

A  direct  extension  of  this  work  would  include  damping  effects  (colli- 
sional  or  Landau).  Inclusion  of  a  phenomenological  damping  term  (damping  rate 
re)  in  the  fluid  equations  for  the  electrons  simply  replaces  ui  by  (<d  +  iTe)  so 
Hill's  method  of  solving  the  infinite  determinantal  dispersion  equation  is 
still  applicable.  Preliminary  calculations  show  no  qualitative  changes,  aside 
from  the  expected  diminution  of  y.  When  Te  exceeds  the  ion  acoustic  frequency 
the  earlier  calculation1  shows  that  the  threshold  with  two  pumps 
can  be  considerably  lower  than  with  a  single  pump  and  it  would  be  interesting 
to  explore  the  behavior  of  the  growth  rates  above  threshold  in  that  case.  An 
analysis  similar  to  that  carried  out  here  could  also  be  used  for  any  of  the 
many  parametric  instabilities  associated  with  magnetized  plasmas.  For  appli¬ 
cations  to  ionospheric  plasmas,  it  would  be  important  to  include  the  effects 
of  density  gradients;  for  example,  each  pump  may  give  rise  to  its  own  decay 
instabilities  albeit  at  different  altitudes,  and  these  would  interact  with  the 
Mathieu  type  of  instabilities  examined  here.  Finally,  a  nonlinear  treatment 
of  two  pump  excitation  would  give  a  more  realistic  prediction  of  the  actual 
line  shapes. 
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8. 


Stability  regions  in  the  plane  of  the  variables  A  anil  B  defined  by 
(41).  The  four  roots  of  the  biquadratic  equation  (40)  for  y  =  ui  -  v> 
have  the  following  characteristics  in  the  respective  regions  of  the  /2 
A-B  plane:  l)  two  real  roots  and  two  conjugate,  purely  imaginary 

roots,  ll)  two  pairs  of  conjugate,  purely  imaginary  roots.  III)  two 
pairs  of  complex  conjugate  roots,  IV)  two  pairs  of  equal  and  oppo¬ 
site  real  roots. 


Stability  regions  in  parameter  space  as  obtained  from  tbe  approxi¬ 
mate  solution  (42).  The  shaded  regions  indicate  the  unstable  zones 
and  the  Roman  numerals  correspond  to  the  labelling  in  Fig.  2.  The 
dotted  line  In  Fig.  Id  is  the  trajectory  In  this  stability  plane 


to  the  single  pump  growth  rate  for  vj  »  1.1  and  gj  =  .346 


ion  acoustic  and  punp  modes  as  a  function  of  vj  with  gj  ■  .346  and 
and  various  values  of  g2.  V£  as  noted. 


Ion  density  n^t)  for  parameters  corresponding  to  the  nulls  of 
Fig.  4  with  gi  =  g2  =  .346,  vL  =  1.1. 

(a)  corresponds  to  null  1  in  Fig.  4  (v2  =  -Vl);  (b)  corresponds  to 
null  2  in  Fig.  4  (vj  -  v2  =  2);  (c)  corresponds  to  null  3  in  Fig.  , 

/.i.  ...  ° 


2 

Fig  •  12:  Electron  density  spectrum  jne(<j)|  for  gj  »  g 2  *  *346,  V1  “  l»l 

plotted  as  a  function  of  (id  -  for  various  choices  of  V2*  ^ 

values  of  (id  -  1 d^)/^  corresponding  to  the  two  pumps  are  indicated 
by  the  dotted  lines. 
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